Infinite series of magnetization plateaus in the frustrated dimer-plaquette chain 
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The dimer-plaquette chain undergoes a first order transition driven by frustration to an orthogonal 
dimer ground state. We analyze the magnetization curve in this orthogonal dimer regime. Besides of 
the recently found magnetization plateaus at m = 1/4 and the m = 1/2 we find an infinite sequence 
of plateaus between m = 1/4 and m = 1/2. These additional plateaus belong to magnetizations 
m{n) = n/{2n + 2), where n is a positive integer number. 



The frustrated dimer-plaquette chain (or orthogonal 
dimer chain) introduced in pi has attracted attention 
because of its first order quantum phase transition to a 
product dimer ground state and its close relation to the 
two-dimensional spin model for SrCu2 (603)2 ^. Moti- 
vated by the recent discovery of magnetization plateaus 
in SrCu2 (603)2 Koga et. al. ||] investigated the mag- 
netization curve for the one-dimensional model. They 
found magnetization plateaus at m ~ 1/4 and m ~ 1/2 
of the full moment near the phase transition point us- 
ing exact diagonalization and DMRG methods. In order 
to reproduce their results we studied the magnetization 
curve using large-size exact diagonalization (up to N=40 
sites) and analytical considerations based on the product 
nature of a certain class of eigenstates of the orthogonal 
dimer chain [n| . We found clear evidence of further mag- 
netization plateaus between m = 1/4 and m — 1/2. In 
what follows we will explain our results for parameters 
J = 1 and j = J' I J = 0.7 (for notations see Fig. 1 and 
Refs. 0,^), i.e., in the zero- field regime the ground state 
is a dimer product state, but the argumentation is more 
general and can be applied in principle to the (b)-phasc 
of the phase diagram presented in Fig. 9 of Q. 

To show that we have not only a single jump between 
the m = 1/4 and m = 1/2 plateaus we need two facts. 

First we have to know the lowest energy of eigenstates 
with m = 1/4 and m = 1/2. The energy for m = 1/4 is 
exactly known as AE/N = {Ei - 3/4)/2 = -1.25595 (see 
table H) , because the corresponding eigenstate is a simple 
product state. For m = 1/2 the ground-state energy can 
be estimated very accurately to iE/N = —0.68059(1) by 
finite-size extrapolation. The table || shows exact diag- 
onalization (ED) results for E{m = 1/2) up to A^ = 40 
sites. Obviously, the energy per site increases with in- 
creasing chain length for an even number of plaquettes, 
whereas the energy per site decreases with increasing 
chain length for an odd number of plaquettes. Both 
energies meet each other at A^ ^ 00. We note, that 
for the considered maximal odd and even chain lengths 
of A^ = 36 and A^ = 40 already the first five digits of 
E{m = 1/2) coincide. 

Second we have to find at least one state with mag- 
netization m' between m — 1/4 and m = 1/2 having an 
energy below the straight line connecting the m — 1/4 
and m = 1/2 points in the E{m) diagram. Then with 



increasing field the system would occupy next after the 
state with to = 1/4 the state with m' and not that one 
with TO, = 1/2 > TO,'. According to [|l| low-lying states 
of certain m belong to the class of product eigenstates, 
for which in every (fc -I- l)-th plaquette the spins ^i and 
53 form a vertical dimer singlet [Si + 6*3)^ ~ (i.e., 
the quantum number 6*13 in Koga's notation is zero). 
The finite strips between two vertical dimer singlets con- 
tain k vertical dimer triplets (6*13 = 1). Since these 
strips are decoupled from each other we call them frag- 
ments of length k here where such a fragment consist of 
Nk = 4fc + 2 spins. The energies E of those eigenstates 
can be expressed by the energies of the finite fragments 
Ek which can be calculated exactly up to a maximum 
fragment length of fc = 9 (see table ||). For instance, the 
energy E of states consisting of N/{Ak + 4) fragments of 
identical length k separated by one vertical dimer singlet 
each is 4^/A^ = {Ek - 3/4)/(fc + 1). 

From the exact numerical data of table H it can be 
found that all of these states have energies below the 
straight line connecting the m = l/A and tti = 1/2 points 
in the E{m) diagram. That clearly shows, that we have 
not only a single jump between the m — 1/4 and the 
m = 1/2 plateau. 

Let us now present a more detailed analysis of the mag- 
netization process between to, = 1/4 and m = 1/2. We 
start from the to, = 1/4 eigenstate, which is a product 
state of identical fragments of length k — 1 each hav- 
ing total spin Sk=i = 1 leading to a total spin of the 
chain Schain = Sk ■ N/{Ak -|- 4) = N/8. To increase the 
magnetization we have to insert one fragment of length 
k = 2 with total spin 5*^=2 = 2. The next eigenstates 
with higher magnetizations can be generated by replac- 
ing further fragments of length fc = 1 by fragments of 
length k — 2. Due to the product nature of these eigen- 
states the energy necessary to increase the magnetization 
is constant and defines the critical field for leaving the 
TO = 1/4 plateau. After all fragments of length fc = 1 
are replaced by fragments of length fc = 2 the system has 
reached a state with Schain = Sk • A^/(4fc + 4) — N/6, i.e., 
a TO = 1/3 plateau. The energy of this to, = 1/3 eigen- 
state is 4S/Af = (£;2-3/4)/3 = -1.076056 (see table Q). 
This fragmented structure of the relevant eigenstates is 
also responsible for the finite-size effects, since for finite 



chains the sphtting of the chain in fragments must be 
compatible with the total length of the chain. To get 
higher magnetizations all fragments of length k — 2 and 
Sk=2 = 2 must now be replaced successively by fragments 
of length fc = 3 and 5^=3 = 3 and so on. This scenario 
generates a series of plateaus at m = k/{2k + 2) where 
the width of the plateaus decreases rapidly, because the 
energy difference between fragments of length k and fc + 1 
is decreasing with increasing k. The critical field between 
two plateaus can be calculated from the energies of the 
fragments Ei^ 



hc{k) ^ kEk - (fc + l)Ek-i + 3/4. 



(1) 



The energies of the fragments Ek can be found in table 
p. Figure gives an overview about the resulting mag- 
netization curve. In contrast to the corresponding Fig .4 
in I3I we get a staircase like infinite series of plateaus. 
We emphasize, according to the linear programming the- 
orem (see, [H and references therein), fragmented states 
with nonuniform fragment lengths have higher energies 
than corresponding states built by fragments of identi- 
cal length. Furthermore our numerical data clearly show 
that the energy of the above considered fragmented states 
is significantly below the energies of competing low-lying 
non-fragmented states (i.e. states where all quantum 
numbers 5*13 = 1). 

We notice, that our result is in accordance with the 
general rule of Oshikawa et al. ^ that n(s — m) has to 
be an integer. In our case the period of the ground state 
n is 4(fc-|- 1), the spin s is one half and the magnetization 
per spin fh is fc/(4fc -I- 4) (the magnetization per site fh 
corresponds to on half of the quantity m used in this 
comment). 

Finally we mention, that the magnetic phase diagram 
on the (j-h) plane presented in Fig. 9 in ||^ has to be 
extended taking into account the additional plateaus at 
TO = fc/(2fc + 2). 

In conclusion, we have found an infinite series of mag- 
netization plateaus in a frustrated one-dimensional quan- 
tum spin system which is closely related to the possibil- 
ity of fragmentation of the considered chain. Though the 
possibility of fragmentation has been observed also for 
frustrated two-leg spin ladder PM], in the spin ladder 
system only very simple ground states seem to be rel- 
evant for its magnetization curve. Hence the sequence 
of plateaus with rational m{k) = fc/(2fc -I- 2) found in 
this paper seems to be a novel property of quantum spin 
systems, not found so far. 
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FIG. 1. The orthogonal dimer chain. 
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FIG. 2. Magnetization m of the orthogonal-dimer spin 
chain for j — 0.7 for the thermodynamical limit (A*' = co). 
The curve m > 1/2 is only a guide for the eyes. The dashed 
line is the exact diagonalization result of the finite system 
AT = 32. 



TABLE I. Finite-size dependence of the energy of m = 1/2 
states E{m = 1/2) for j = 0.7 up to Af = 4Np = 40 sites. 
Np is the number of plaquettes (vertical dimers) in the chain 
with periodic boundary conditions. 
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Np 


E/Np 


Np 


E/Np 


1 


-0.50000000 


2 


-0.70313092 


3 


-0.67602844 


4 


-0.68187420 


5 


-0.68023225 


6 


-0.68071087 


7 


-0.68055572 


8 


-0.68060683 


9 


-0.68058916 


10 


-0.68059535 



TABLE II. Energies of fragments of the orthogonal-dimer 
chain of length k and total spin Sk = k for j = 0.7 and up to 

ATfc = 4fc + 2 = 38 sites. 



k 


Ek 


k 


Ek 





-0.7500000 


5 


-4.5289991 


1 


-1.7619002 


6 


-5.2096756 


2 


-2.4781677 


7 


-5.8902882 


3 


-3.1658665 


8 


-6.5708863 


4 


-3.8480439 


9 


-7.2514810 



